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1. Introduction 

An existing dialogue between number theory and dynamical systems 
is advanced. A combinatorial device gives necessary and sufficient con- 
ditions for a sequence of non-negative integers to count the periodic 
points in a dynamical system. This is applied to study linear recurrence 
sequences which count periodic points. Instances where the p-parts of 
an integer sequence themselves count periodic points are studied. The 
Mersenne sequence provides one example, and the denominators of the 
Bernoulli numbers provide another. The methods give a dynamical 
interpretation of many classical congruences such as Euler-Fermat for 
matrices, and suggest the same for the classical Kummer congruences 
satisfied by the Bernoulli numbers. 

Let Mn = 2" — l,n > 1 denote the n-th term of the Mersenne se- 
quence {Mn). This sequence is of interest in number theory because it 
is expected to contain infinitely many prime terms, and in dynamics 
because it counts the periodic points in the simplest expanding dynam- 
ical system. Let T : — ^ be the squaring map T{z) = z"^, and let 
Per„(T) denote the set of points of period n under T, that is the set 
of solutions of the equation T^{z) = z. Then it is easy to check that 
|Per„(T)| =M„. 

Other classical sequences arise in a similar way. Let denote the n- 
th term of the Lucas sequence 1, 3, 4, 7 ... , and let X denote the set of 
all doubly- infinite strings of O's and I's in which every zero is followed 
by a 1, and T : X ^ X the left shift defined by {Tx)n = Xn+i- Then 
\Vein{T)\ = Ln. 

The Lehmer-Pierce sequences (generalising the Mersenne sequence; 
see IQ) also arise in counting periodic points. Let f{x) denote a monic. 
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integral polynomial with degree d > 1 and roots ai, . . . , a^- Define 

An(/) = n 



la" 



which is non-zero for n > 1 if no is a root of unity. When f{x) = 
X — 2, we obtain A„(/) = M„. Sequences of the form (A„(/)) were 
studied by Pierce and Lehmer with a view to understanding the spe- 
cial form of their factors, in the hope of using them to produce large 
primes. In dynamics they arise as sequences of periodic points for 
toral endomorphisms. Let X = T'^ denote the (i-dimensional additive 
torus. Then the companion matrix Af of f acts on X by multiplication 
mod 1, T{x) = AfX mod 1. It requires a little thought to check that 
I Per„(T)| = A„(/) under the same ergodicity condition that no is 
a root of unity (see 0). Notice that the Lehmer-Pierce sequences are 
the absolute values of integer sequences which could have mixed signs. 

Our final examples illuminate the same issue of signed sequences 
whose absolute value counts periodic points. The Jacobsthal-Lucas 
sequence Rn = |(— 2)" — 1| counts points of period n for the map 
z I— >• z~'^ on S^. The sequence Sn = |2" -|- (— 3)"| counts periodic points 
in a certain continuous automorphism of a 1-dimensional solenoid, see 
m or 



Following [T^, call a sequence of non-negative integers realisable 
if there is a set X and a map T : X ^ X such that m„ = | Per„(T)|. 
The examples above were of sequences that are realisable by continuous 
maps of compact spaces; it turns out that any realisable sequence is in 
fact realisable by such a map. 

It is natural to ask what is required of a sequence in order that it 
be realisable. For example, could the Fibonacci sequence, the more 
illustrious cousin of the Lucas sequence, be reahsed in this way? The 
answer is no, and a simple proof will follow in the next section. In fact 
a sequence of non-negative integers satisfying the Fibonacci recurrence 
is realisable if and only if it is a non-negative integer multiple of the 
Lucas sequence (see |jT2| and ||13|)- 



The statements of the main theorems now follow. For the first, note 
that if (un) is any sequence of integers, then it is reasonable to ask if 
the sequence {\un\) of absolute values is reahsable. For example, the 
sequence 1, — 3, 4, — 7, . . . is a signed linear recurrence sequence whose 
absolute values are realisable. 

The first theorem gives a generalisation of the observation about 
realisable sequences which satisfy a linear recurrence relation. The 
definitions are standard but they will be recalled in the next section. 
Recall that the C-space of all solutions of a binary recurrence relation 
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has dimension 2. The realisable subspace is the subspace generated by 
the reahsable solutions. For the Fibonacci recurrence, the reahsable 
subspace has dimension 1 and is spanned by the Lucas sequence. 

Theorem 1.1. Let u„ denote the n-th term of an integer sequence 
which satisfies a non- degenerate binary recurrence relation. Let A de- 
note the discriminant of the characteristic polynomial associated to the 
recurrence relation. Then the realisable subspace has 

1. dimension if A < 0, 

2. dimension 1 if A = or A > and non-square, and 

3. dimension 2 if A > is a square. 

Theorem |1.1| surely has a generalisation to higher degree which char- 
acterises the realisable subspace in terms of the factorisation of the 
characteristic polynomial /. The second theorem is a partial result in 
that direction, giving a restriction on the dimension of the realisable 
subspace under the assumption that the characteristic polynomial has 
a dominant root. 

Theorem 1.2. Let f denote the characteristic polynomial of a non- 
degenerate linear recurrence sequence with integer coefficients. If f is 
separable and has I irreducible factors and a dominant root then the 
dimension of the realisable subspace is < I. 

It is not clear if there is an exact result, but the deep result of 
Kim, Ormes and Roush ^ on the Spectral Conjecture of Boyle and 
Handelman gives a checkable criterion for a given linear recurrence 
sequence to be realised by an irreducible subshift of finite type. 

Example 1.3. Consider the sequences which satisfy the Tribonacci 
relation 

The sequence 1,3,7,11,21,... satisfies (0) and is realisable. This 
is the sequence of traces T(A"), where A is the companion matrix to 
f{x) = — — X — 1., 

A = 




For an explanation of this remark, turn to the proof of Corollary ^A . 
Theorem [1.2| says that any realisable sequence which satisfies (0) is a 
multiple of this one. 

The third theorem consists of a pair of examples. Given a sequence 
Un and a prime p, write [un]p for the p-part of We say a sequence 
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is locally realisable at p if the sequence [un]p is itself realisable. We say 
the sequence is everywhere locally realisable if it is locally realisable at 
p for all primes p. If a sequence is everywhere locally realisable, then 
for each n > 1, [un\p = 1 for all but finitely many p, and it is realisable 
by Corollary |2.2| . We will sometimes use the term globally realisable 
for a sequence when we wish to emphasize the distinction with local 
realisability. Consider the Bernoulli numbers, which are defined by the 
formula 

t °° t" 

n=0 

Then _B„ G Q for all n, and -B„ = for all odd n > I. 

Theorem 1.4. 1. Any Lehmer-Pierce sequence is everywhere locally, 
and hence globally, realisable. 
2. Let bn denote the denominator of for n > 1. Then is 
everywhere locally, and hence globally, realisable. 

The maps in Theorem p. .41 are endomorphisms of groups. Theorem 



ni and Lemma |2.1| suggest a dynamical interpretation of composite 



versions of the classical Kummer congruences; see section ^ below. 



2. Combinatorial dynamics 

We begin with a simple remark that shows the Fibonacci sequence is 
not realisable. No map can have 1 fixed point and exactly 2 points of 
period 3, as any point of least period 3 must have an orbit of length 3 
comprising points all of period 3. More generally, for any prime p, the 
number of non-fixed points of period p must be divisible by p because 
their orbits occur in cycles of length p. Using a generalisation of this 
kind of reasoning, the following characterisation emerges. 

Lemma 2.1. Let u = (m„) be a sequence of non-negative integers, 
and let u * n denote the Dirichlet convolution of u with the Mobius 
function ft. Then u is realisable if and only if {u * = mod n 

and (u * fi){n) > for all n > 1. 

To see why this holds, notice that the set of points of period n is 
the disjoint union of the set of points of least period d for d running 
through the divisors of n, and the number of points with least period d 
is a mutliple of d. The Dirichlet convolution is the usual definition from 
analytic number theory: = 1, fj,{n) = unless n is square-free and 
fj,{n) = (—1)'^ if n is the product of k distinct primes, and the Dirichlet 
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convolution of an arithmetical function g with fi is given by 
{g * fi){n) = ^n{d)g{n/d). 

d\n 

Finally, the result is obtained using the Mobius inversion formula. For 
brevity, write u* = (n * fi){n) for n > 1. 

Corollary 2.2. The sum and product of two realisable sequences are 
both realisable. 

This may be seen either using elementary properties of the Dirichlet 
convolution or using the realising maps: if u and v are realisable, then 
the Cartesian product of the realizing maps realises (unVn), while the 
disjoint union realises (u„ + Vn)- 



Notice that if n = j?*", for a prime p and r > an integer. Lemma |2?1 
requires that 

Upr = Upr-1 mod p^ 
for any realisable sequence u. 

Corollary 2.3. Let a denote a positive integer and let p and r be as 
above. Then 

= mod p^' . 

Proof. This is the statement of the Euler-Fermat Theorem, which may 
be seen because the sequence u„ = a"' is realisable. For example, the 
left shift T on {0, 1, . . . , a - 1}^ has | Per„(T) | = a". □ 

This kind of observation — that periodic points in full shifts give 
simple proofs of many elementary congruences — is folklore; indeed 
the paper gives a rather complicated proof of Euler-Fermat using a 
dynamical system. 

Lemma |2.1| does more with no additional effort. The following is a 
generalisation of the Euler-Fermat Theorem for integral matrices which 



will be used in the proof of Theorem |LT 



Corollary 2.4. Let A denote a square matrix with integer entries and 
let p and r be as above. Then 

trace(y4^'^) = trace (A^*^ ^) mod p^. 

Proof. It is sufficient to assume A has non-negative entries, since any 
matrix has such a representative mod p'". For non-negative entries, 
(trace(A")) is realisable: Let Ga be the labeled graph with adjacency 
matrix A and Ta the edge-shift on the set of labels of infinite paths on 
Ga- Then the number of points of period n for this system is trace(yl") 
(see |11|] for the details). □ 
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We now state the consequences of Lemma ^]T| in their most general 
form for matrix traces. 

Corollary 2.5. Let A denote a square matrix with integer entries and 
let An denote the sequence trace(y4'^). Then for all n > 1 

A* = mod n. 



Before the proof of Theorem [L^, we begin with some notation (for a 



lively account of the general properties of linear recurrence sequences, 
see 10). Suppose we are given a binary recurrence sequence u = (u„). 
This means that Ui and U2 are given as initial values with subsequent 
terms defined by a recurrence relation 

(2) Un+2 = BUn+1 - CUn- 

The polynomial f{x) = — Bx + C is the characteristic polynomial 
of the recurrence relation. We will write 

1 



- y-c B 

for the companion matrix of /. The zeros ai and a2 of /, are the 
characteristic roots of the recurrence relation. The assumption on non- 
degeneracy means that ai/a2 is not a root of unity. The discriminant 
of the recurrence relation is A = i?^ — AC. Of course, if A = then 
the roots of / coincide, if A < the roots are non-real and distinct, if 
A > is a square then the roots are rational and in the other case, the 
roots are real and distinct but irrational. 

The general solution of the recurrence relation in these cases is as 
follows: 

A = 0: Un = (71 + 72'^)«i (here ai = 02). 
A 7^ 0: Un = 7i"? + 72^2 • 



Proof, (of Theorem |LT 



Assume first that A = 0, and let p denote any prime which does not 
divide 72. Then the congruence in Corollary |2.4| is plainly violated at 
n = p unless 72 = 0. In that case, |7iQ;i | is realisable and the space 
this generates is 1-dimensional. 

If A > is a square, then the roots are rationals and plainly, must 
be integers. We claim that for any integers 71 and 72, the sequence 
|7ia;" -|- 72^2 I is realisable. In fact (up to multiplying and adding full 
shifts) this sequence counts the periodic points for an automorphism 
on a one-dimensional solenoid, see Q or |T^ . 



The two cases where A 7^ is not a square are similar. Write 
a = e + /a/a for one of the roots of / and let K = Q{a) denote 
the quadratic number field generated by a. Write T^iq : K Q for 
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the usual field trace. The general integral solution to the recurrence 
is Un = Ti^|(Q((a + b^/A)a"'), where a and b are both integers or both 
half-odd integers. Write t>„ = Ti^|Q(aa") and Wn = TK\Q{b\/Aa^). 
Now Vn = atrace(ylj), where Af denotes the companion matrix of /. 
Hence it satisfies Vp = Vi mod p for all primes p by Corollary p.4| . 

Let p denote any inert prime for K. The residue field is isomorphic 
to the field Fp2. Moreover, the non-trivial field isomorphism restricts 
to the Frobenius at the finite field level. Reducing mod p gives the 
congruence 



Wp 



T (yAa^ = \^aP - ^Aa mod p. 



Thus, Wp = —Wi mod p for all inert primes p. On the other hand, 
Vp = f 1 mod p for all inert primes p. 

If \un\ is realisable then \up\ = \ui\ mod p by Corollary |2]^. If 
Up = —Ml mod p for infinitely many primes p then Vp + Wp = Vi — Wi = 
—vi — wi mod p. We deduce that p\vi for infinitely primes and hence 
vi = 2ae = 0. We cannot have e = by the non-degeneracy, so a = 0. 
If Up = Ui mod p then, by a similar argument, we deduce that bf = 0. 
We cannot have / = again, by the non- degeneracy so 6 = 0. This 
proves that when A 7^ is not a square, the realisable subspace must 
have rank less than 2. 

Suppose firstly that A > 0. We will prove that the rank is precisely 
1. In this case, there is a dominant root. If this root is positive then 
all the terms of Un are positive. If the dominant term is negative 
then the sequence of absolute values agrees with the sequence obtained 
by replacing a by —a and the dominant root is now positive. In the 
recurrence relation (^) C = Nx|Q(a), the field norm, and B = T;^|Q(a). 
We are assuming B > 0. If C < then the sequence Un = T{A^) is 
realisable because the matrix Af has non-negative entries. If C > 
then we may conjugate ^/ to such a matrix (this leaves the sequence 
of traces invariant). To see this, let E denote the matrix 

1 0^ 
k 1 



E 
Then 



E~'^fE- [Bk-%-C B-k 

If B is even, take k = B/2. Then the lower entries in E~^AfE are 
- 4C)/4 = A/4 > and B/2 > 0. li B is odd, take k = {B + l)/2. 
Then the lower entries are (5^ - 1 - 4C)/4 = (A - l)/4 > and 
{B — l)/2 > 0. In both cases we have conjugated Aj to a matrix with 
non-negative entries. Since we know that the sequence of traces of a 
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matrix with non-negative entries is realisable, we have completed this 
part of the proof. 

Finally, we must show that when A < 0, both sequences Vn and 
Wn are not realisable in absolute value. Assume a 7^ 0, and then 
note that Vi = 2ae 7^ by the non-degeneracy assumption. For all 
primes p we have Vp = Vi by the remark above. Since the roots ai 
and a2 are complex conjugates, |q;i| = \a2\. Let /3 = 2^ arg(a;i/a2); 
f3 is irrational by the non-degeneracy assumption. The sequence of 
fractional parts of pf3, with p running through the primes, is dense 
in (0,1) (this was proved by Vinogradov |T6[; see |T^ for a modern 



treatment). It follows that there are infinitely many primes p for which 
VpVi < 0. Therefore, if is realisable then it satisfies Vp = Vi mod 
p and —Vp = Vi mod p for infinitely many primes. We deduce that 
f 1 = which is a contradiction. With w„ we may argue in a similar 
way to obtain a contradiction to t^i 7^ 0. If is realisable then 
Lemma P?T] says \wp2\ = \wp\ = \wi \ for all primes p. Arguing as before, 
Wp2 = wi for both split and inert primes. However, the sequence {p'^P}, 
p running over the primes, is dense in (0,1). (Again, this is due to 



Vinogradov in [|l^] or see for a modern treatment. The general case 
of {F{p)}, where F is a polynomial can be found in [0.) We deduce 
that Wp2Wi < for infinitely many primes. This means Wp2 = Wi 
mod p and Wp2 = —wi mod p infinitely often. This forces ti^i = - a 
contradiction. □ 



3. Proof of Theorem 1.2 



The proof of Theorem |r]2| uses the methods introduced in the proof 
of Theorem |1.1| . 

Proof. Let d denote the degree of /. In the first place we assume 
1 = 1, thus / is irreducible. The irreducibility of / implies that the 
rational solutions of the recurrence are given by m„ = T Kiqi'^d"') , where 
K = Q(a), and 7 G if. We write ■yi,ai,i = 1, . . . ,d for the algebraic 
conjugates of 7 and a. The dominant root hypothesis says, after re- 
labelling, |q;i| > \ai\ for z = 2,...,d. We will show that if u„ is 
realisable then 7 G Q. 

Let p denote any inert prime. If p is sufficiently large, the dominant 
root hypothesis guarantees that Up, . . . , Upd will all have the same sign. 



Using Lemma |2.1| several times, we deduce that 

Up = Up2 = ■ ■ ■ = Upd = ±ui mod p. 
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Therefore Wp + ■ ■ ■ + Upd = ±dui mod p, the sign depending upon the 
sign of ui. However, 

Up-\ hUpd = T;^|q(7) Ti^|Q(") mod p. 

We deduce a fundamental congruence 

Tk|q(7) Tx|Q(a) = ±dTK\Q{ja) mod p. 

Since this holds for infinitely many primes p, the congruence is actually 
an equality, 

(3) T;^|q(7) Tx|Q(a) = ±dTK\Qi'ja). 

The next step comes with the observation that if m„ is realisable then 
Urn is realisable for every r > 1. Thus equation @ now reads 

(4) T;^|q(7) T/<|Q(a'") = ±rfTx|Q(7«'')- 

Dividing equation (^) by a[ and letting r ^ oo we obtain the equation 

Ti^|Q(7) = ±c?7i- 

This means that one conjugate of 7 is rational and hence 7 is rational. 

The end of the proof in the case / = 1 can be re-worked in a way 
that makes it more amenable to generalisation. The trace is a Q-linear 
map on K so its kernel has rank d — 1. Thus every element 7 of 
can be written g + 70 where g G Q and Tx|(q(7o) = 0. Noting that 
Tx|(q((?) = dq and cancelling d, this simply means equation (Q) can be 
written 

Ur = ±qTK\Q{a''), 

for all r > 1 confirming that the realisable subspace has rank < 1. 

The general case is similar. Each of the irreducible factors of / 
generates a number field Kj,j = of degree dj = [Kj : Q]. The 

solutions of the recurrence look like 

I 

where each G Kj. Let L denote the compositum of the Kj. Using 
the inert primes of L and noting that each is inert in each Kj, we 
deduce an equation 

I ^ I 

(5) J. ^K.nilj) Tx,|Q(ai) = ±d'Yi:K,\Q{ljaj). 
j=i J j=i 
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As before, replace a,- by a^, and cancel d so that 



1 



j=l J 



Each 7j can be written 7j = qj + joj, where TRjiqiloj) = 0. Noting 
that Tx-iQigj) = djQj we deduce that 



Ur 



which proves that the realisable subspace has rank < I. □ 
4. Proof of Theorem O 



It is sufficient to construct local maps Tp : Xp — >■ Xp for each prime 
p. Then Corollary |2.2| guarantees a global realisation by defining 

T = l[Tp on X = l[Xp 

p p 

If the maps Tp are group endomorphisms then the map T is a group 
endomorphism. 

Proof. As motivation, consider the Mersenne sequence. For each prime 
p, let Up C denote the group of all pth power roots of unity. Define 
the local endomorphism Sp : x on Up. Then | Per„(S'p) | = [2*^ — l]p 

so Sp gives a local realisation of the Mersenne sequence. 

An alternative approach is to use the S'-integer dynamical systems 
from 0: for each prime p, define Tp to be the automorphism dual to 
X i-H> 2x on Z(p) (the localisation at p). Then by 0, 

I Per„ (Tp) I = H |2"-1|,= [2"-l]p 

by the product formula. This approach gives a convenient proof of the 
general case. We may assume that the polynomial / is irreducible; let 
K = Q(^) for some zero of /. Then for each prime p, let S comprise all 
places of K except those lying above p, and let Tp be the S'-integer map 
dual to X ^ ^x on the ring of S'-integers in K. Then by the product 
formula 

|Per„(Tp)| = (nir-l|.)"' = [A„(/)], 

v\p 

as required. 

For the Bernoulli denominators, define Xp = Fp = Z/pZ. For p = 2 
define Tp to be the identity. For p > 2, let ^fp denote an element 
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of (multiplicative) order [p — l)/2. Define Tp : Xp Xp to be the 
endomorphism Tp(x) = QpX mod p. Plainly |Per„(Tp)| = p if and only 
if p — l|2n; for all other n, |Per„(Tp)| = 1. The von Staudt-Clausen 
Theorem (0, 0) states that 

B2n + V - G Z, 

^-^ p 

where the sum ranges over primes p for which p—1 1 2?T,. Thus |Per„(Tp)| = 
max{l, |-B2n|p} and this shows the local realisability of the Bernoulli 
denominators. □ 

The following statements seem plausible upon numerical investiga- 
tion. 

1. the denominators 

12,120,252,240,132,32760,... 

of -B2n/2^ form a sequence that is everywhere locally realisable; 

2. the numerators 

1,1,1,1,1,691,1,3617,43867,... 

of i?2n/2n form a realisable sequence that is not locally realisable 
at the irregular primes 37, 59, 67, 101, 103, 131, 149, 157, .... 

3. the denominators 

24, 240, 504, 480, 264, 65520, 24, 16320, . . . 

of B2n/ 4n form a realisable sequence that is not locally realisable 
at the primes 2, 3, 5, 7, 11, 13 but seems to be locally realisable for 
large primes. 

Taking these remarks together with n = p'' in Lemma |2.1| , suggests 
a dynamical interpretation of the Kummer congruences. These are 
stated now, for a proof see |0. 

Theorem 4.1. If p denotes a prime and p — 1 does not divide n then 
n = n' mod {p — l)p^ implies 

/I n-l\Bn n'—l\Bn' j r+1 

[I — p ) — = [1 — p ) mod p . 

n n' 
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